The Scalar Curvature of a Causal Set 
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A one parameter family of retarded linear operators on scalar fields on causal sets is introduced. 
When the causal set is well approximated by 4 dimensional Minkowski spacetime, the operators are 
Lorentz invariant but nonlocal, are parametrised by the scale of the nonlocality and approximate 
the continuum scalar D'Alembertian □ when acting on fields that vary slowly on the nonlocality 
scale. The same operators can be applied to scalar fields on causal sets which are well approximated 
by curved spacetimes in which case they approximate □ — ^i? where R is the Ricci scalar curvature. 
This can used to define an approximately local action functional for causal sets. 

PACS numbers: 04.60.Nc,02.40.-k,11.30.Cp 
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The coexistence of Lorentz symmetry and fundamen- 
tal, Planck scale spacetime discreteness has its price: one 
must give up locality. Since, if our spacetime is granu- 
lar at the Planck scale, the "atoms of spacetime" that 
are nearest neighbours to a given atom will be of or- 
der one Planck unit of proper time away from it. The 
locus of such points in the approximating continuum 
Minkowski spacetime is a hyperboloid of infinite spa- 
tial volume on which Lorentz transformations act tran- 
sitively. The nearest neighbours will, loosely, comprise 
this hyperboloid and so there will be an infinite number 
of them. Where curvature limits Lorentz symmetry, it 
may render the number of nearest neighbours finite but 
it will still be huge so long as the radius of curvature is 
large compared to the Planck length. Causal set theory 
is a discrete approach to quantum gravity which embod- 
ies Lorentz symmetry y] and exhibits nonlocality of 
exactly this form 0, Q . 

Nonlocality looks to be simultaneously a blessing and 
a curse in tackling the twin challenges that any funda- 
mentally discrete approach to the problem of quantum 
gravity must face. These are to explain (1) how the fun- 
damental dynamics picks out a discrete structure that 
is well approximated by a Lorentzian manifold and (2) 
why, in that case, the geometry should be a solution of 
the Einstein equations. This is often referred to as the 
problem of the continuum limit but in the context of a 
fundamentally discrete theory in which the discreteness 
scale is fixed and is not taken to zero but rather the ob- 
servation scale is large, it is more accurately described 
as the problem of the continuum approximation. 

Consider first the problem of recovering a continuum 
from a quantum theory of discrete manifolds. (We adopt 
this term following Riemann Q and use it to refer to 
causal sets, simplicial complexes, graphs, or whatever dis- 
crete entities the underlying theory is based on.) When- 
ever a background principle or structure in a physical the- 
ory is abandoned in order to seek a dynamical explana- 
tion for that structure, the state we actually observe be- 
comes a very special one amongst the myriad possibilities 
that then arise. The continuum is just such a background 
assumption. In giving it up, generally one introduces a 



space of discrete manifolds in which the vast majority 
have no continuum approximation. There will therefore 
be a competition between the entropic pull of the huge 
number of noncontinuum configurations - choose one uni- 
formly at random and it will not look anything like our 
spacetime - and the dynamical law which must suppress 
the contributions of these nonphysical configurations to 
the path integral. The following general argument shows 
that a local dynamics for quantum gravity will struggle 
to provide the required suppression. Consider the parti- 
tion function as a sum over histories in which the weight 
of each discrete manifold is e~'^ where S is the real Wick 
rotated action. As we increase the observation scale, the 
sum will be over discrete manifolds with an increasing 
number, A'^, of atoms. If the action is local - which in 
a discrete setting translates to it being a sum over con- 
tributions from each atom - then it will grow no faster 
than N times some constant, a, and so each weight is 
no smaller than e~"^. If the number of discrete mani- 
folds with N atoms grows faster than exponentially with 
N, and if the majority of these discrete manifolds are 
not continuumlike then they will overwhelm the parti- 
tion function and the typical configuration will not have 
a continuum approximation. Even when the number of 
discrete manifolds is believed to grow exponentially, en- 
tropy can still trump dynamics as was seen in the lack of 
a continuum limit in the Euclidean dynamical triangula- 
tions programme [6|- |9|. C ausal dynamical triangulations 
do better, see, e.g., |10l - ll3j . by restricting the class of 



triangulations allowed in the sum. 

In the case of causal sets, the number of discrete man- 



ifolds of size N grows as e^ 1^ 14 1 and a local action 
would give causal set theory little chance of recovering 
the continuum. So the nonlocality of causal sets holds out 
hope that the theory has a continuum regime and indeed 
there exist physically motivated, classically stochastic dy- 
namical models for causal sets 15| in which the entropi- 
cally favoured configurations almost surely do not occur 
and those that do exhibit an intriguing hint of manifold- 
like-ness (l6| . 

However, nonlocality poses a danger when it comes to 
the second challenge of recovering Einstein's equations. 
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If we assume that a discrete quantum gravity theory does 
have a 4 dimensional continuum regime, and if the theory 
is local and generally covariant, then the long distance 
physics will be governed by an effective Lagrangian which 
is a derivative expansion in which all diffeomorphism in- 
variant terms are present but higher derivative terms are 
suppressed by the appropriate powers of the Planckian 
discreteness length scale, I: 



ail-^R + a2R^ 



(1) 



where R is the Ricci scalar, ai and 02 are dimensionless 
couplings of order 1, and the dots denote further cur- 
vature squared terms as well as cubic and higher terms. 
The coefficient of the leading term, oq, is also naturally 
of order 1 which would make it 120 orders of magnitude 
larger than its observed value. However, that would also 
produce curvature on Planckian scales and so would not 
be compatible with the assumption of a continuum ap- 
proximation. In a discrete theory, the question of why the 
cosmological constant does not take its natural value is 
the same question as why there is a continuum regime at 
all and we must look to the fundamental dynamics for its 
resolution. Assuming there is a resolution and a contin- 
uum regime exists, locality and general covariance then 
pretty much guarantee Einstein's equations due to the 
natural suppression of the curvature squared and higher 
terms compared to the Einstein-Hilbert term. 

So, Lorentz symmetry and discreteness together imply 
nonlocality, but nonlocality blocks the recovery of general 
relativity, and if causal sets were incorrigibly nonlocal, 
this would be fatal. Suppose, however, that the nonlocal- 
ity were somehow limited to length scales shorter than a 
certain Ik-, which could be much larger than the Planckian 
discreteness scale, but yet have remained experimen- 
tally undetected to date. There is already evidence that 
this is possible and indeed causal sets admit construc- 
tions that are local enough to approximate the scalar 
D'Alembertian operator in 2 dimensional flat spacetime 
[itI [isj . We add to this evidence here by exhibiting a 
family of discrete operators that approximate the scalar 
D'Alembertian in 4 dimensional flat spacetime. Further, 
both the 2D and 4D operators, when applied to scalar 
fields on causal sets which are well described by curved 
spacetimes approximate □ — ^i?, where R is the Ricci 
scalar curvature. We use this to propose an action for a 
causal set which is approximately local. 

We recall that a causal set (or causet) is a locally finite 
partial order, i.e., it is a pair (C, ^) where C is a set and 
^ is a partial order relation on C, which is (i) reflexive: 
X ^ X, (ii) acyclic x^y^x^x — y, and (iii) transitive 
x:^y^z^x^z, for all x,y, z C. Local finiteness is 
the condition that the cardinality of any order interval is 
finite, where the (inclusive) order interval between a pair 
of elements y ^ a; is defined to be I{x, y) ■= {z ^ C\y < 
z ^ x}. We write x ^ y when x < y and x ^ y. We 



call a relation x < y a, link if the order interval I{x, y) 
contains only x and y: they are nearest neighbours. 

Sprinkling is a way of generating a causet from a d- 
dimensional Lorentzian manifold {A4 ,g). It is a Poisson 
process of selecting points in M. with density p so that the 
expected number of points sprinkled in a region of space- 
time volume V is pV . This process generates a causet 
whose elements are the sprinkled points and whose order 
is that induced by the manifold's causal order restricted 
to the sprinkled points. We say that a causet C is well 
approximated by a manifold {AA,g) if it could have been 
generated, with relatively high probability, by sprinkling 
into {M,g). 

We propose the following definition of a discrete 
D'Alembertian, B, on a causet C that is a sprinkling, 
at density p = l^'^, into 4D Minkowski space M**. Let 
: C — ?> R be a real scalar field, then 

B(t>(x) :=^[-0(x) 
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where the sums run over 4 layers Li, i — 1, . . . , 4, 

Li := {y E C : y ^ x and n{x, y) — i ~ 1} (3) 

and n{x,y) := \I{x,y)\ — 2. So, for example, layer Li 
is the set of all elements y that are linked to x and as 
described above, they will be distributed close to a hy- 
perboloid that asymptotes to the past light cone of x and 
is proper time I away from x. This sum will not in gen- 
eral be uniformly convergent if it is over the elements of 
a sprinkling into infinite so we introduce an IR cut- 
off, L >> I, by embedding C in M'^ and summing over 
the finitely many elements sprinkled in the intersection 
of the causal past of x and a ball of radius L centred on 
X. The details of the calculation that shows why 4 layers 
are necessary in 4D will appear elsewhere, however see 
[lil for an explanation of why 3 layers are needed in 2D 
and the conjecture that 4D will require 4 layers. 

Now let (p he a. real test field of compact support on 
M^. If we fix a point x G (which we always take to be 
included in C) and evaluate B(j){x) on a sprinkling into 
M^, its expectation value in this process is given by 

B0(x) :^E(i30(a:)) = ^[-</.(x) 

+ 1 / d'y 0(y) e~« (1 - 9^ -f %e -\e)], (4) 

' JyeJ-{x) 

where ^ l^'^V{x, y), V{x, y) is the volume of the causal 
interval between x and y and there is an implicit cutoff 
L, the size of the support of <j), on the integration range. 

It can be shown that this mean converges, as the 
discreteness scale is sent to zero, to the continuum 
D'Alembertian of (j). 



lim B<j)ix) 
1^0 



n(j){x) 



(5) 
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and that B(j>{x) is well approximated by D(j){x) when 
the characteristic length scale, A, on which 4'{x) varies is 
large compared to L B is therefore effectively sampling 
the value of the field only in a neighbourhood of x of size 
of order I and the mean, at least, of B is about as local 
as it can possibly be, given the discreteness. 

To see roughly how this can happen, notice that the in- 
tegrand in dU is negligible for ^ > a'' where a is such that 

4 

e~" << 1. The significant part of the integration range 
therefore lies between the past light cone of x and the 
hyperboloid ^ = and comprises a part within a neigh- 
bourhood of X of size al - whence the local contribution 
- and the rest which stretches off far down the light cone. 
It is this second part of the range which threatens to in- 
troduce nonlocality but because it can be coordinatized 
by ^ itself and some coordinates 77" on the hyperboloid 
the integration over it will be proportional to 



^4 



(6) 



If is nearly constant over length scale al, the ^ integra- 
tion is close to zero and the contribution is suppressed. 

The fluctuations in B(j){x), however, are a different 
matter: if the physical IR cutoff L is fixed and the dis- 
creteness scale sent to zero, i.e., the number of causet ele- 
ments N grows, simulations show the fluctuations around 
the mean grow rather than die away and B(j){x) will not 
be approximately equal to the continuum \3(j)(x). To 
dampen the fluctuations we follow [3| and introduce an 
intermediate length scale Ik > I and smear out the ex- 
pressions above over this new scale, with the expectation 
that when 1^ » I the inhering averaging will suppress 
the fluctuations via the law of large numbers. Thus we 
seek a discrete operator, Bk, whose mean is given by dU 
but with I replaced by Ik- 



Bk(j){x) 



Veil 



+ i f m e~«(i - 9^ + \e)\ , (7) 

'fc Jye.J-{x) •J 

where now ^ :— l'f,'^V{x, y). Working back, one can show 
that the discrete operator, Bk, with this mean is 



Bk(t){x) 



Veil 

where e = (l/lkY and 
/(n,6) = (l-er 



y-<x 



9en 



, (8) 



1-e (n-2)!(l-e)2 
4e3n! 



3(n-3)!(l-e)3 



(9) 



Bk reduces to B when e — 1. Bk effectively samples 
(h over elements in 4 broad bands with a characteristic 



depth Ik, the bands' contributions being weighted with 
the same set of alternating sign coefficients as in Since 
JT]) is just dH) with I replaced by Ik, the mean of Bk(j){x) is 
close to □0(a;) when the characteristic scale over which 
varies is large compared to Ik- Now, however, numerical 
simulations show that the fluctuations are tamed. Points 
were sprinkled into a fixed causal interval in between 
the origin and t = 1 on the t axis, at varying density p = 
y, where volume V = j^. For each N, 100 sprinklings 
were done and for each sprinkling, Bk4> was calculated at 
the topmost point of the interval for 0=1 and Ik = 0.16. 
For N — 5000, the mean was /i = 9, 35 and the standard 
deviation s.d = 134.8. For N = 10000, fi = -4.00 and 
s.d. = 102.6 and for N = 20000, p = 1.12 and s.d. = 
58.8. These results indicate that the fluctuations do die 
away, as anticipated, as N increases and are consistent 
with the dependence iV^"^/^) ^ Further results will appear 
elsewhere. 

The operators B and Bk derived in both 2D (in [18]) 
and 4D are defined in terms of the order relation on C 
alone and so can be applied to a scalar field on any causet. 
If, therefore, {A4,g) is a (2D or 4D) curved spacetime and 
is a scalar field on Ai, we can compute Bk4>{x) on a 
sprinkling into M and calculate its mean. Let V2 and V4 
be the volumes of the intervals in 2D and 4D respectively, 
^2 := y2{x,y)l~^'^ and ^4 := Vi{x,y)l'i^^ . Then, in the 
presence of curvature. 



B^kU{x) 
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B 



(4), 



{x) = 



1 



(10) 



Vblj. Ik JyeJ-(x) 



(1 - 96 + : 



-Al)^], (11) 



in 2D and 4D respectively. 

These expressions can be evaluated using Riemann 
normal coordinates and in both cases we find 



lim B 



1 



= □ - -R{x, 



(12) 



The limit is a good approximation to the mean when the 
field (f) varies slowly over length scales Ik and the radius 
of curvature r » Ik- 

If the damping of fluctuations found in simulations in 
fiat space are indicative of what happens in curved space 
then, for a fixed large enough IR cutoff, L, the nonlocality 
length scale Ik can be chosen such that I « Ik « L and 
the value of Bk4) for a single sprinkling will be close to the 
mean. If Bk is applied to the constant field </> = — 2, we 
therefore obtain an expression that is close to the scalar 
curvature of the approximating spacetime. 
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In each of 2D and 4D, we can now define a one pa- 
rameter family of candidate actions, S'^p], for a causal 
set, C, by summing B^^—l) over the elements of C, times 
hl'^ to get the units right, times a number of order one 
which in 4D is the ratio of P to 1^, where Ip = VSttGH 
is the rationalized Planck length. When the nonlocality 
length Ik equals the discreteness length I, Bk = B and 
the action, S[C] takes a particularly simple form as an 
alternating sum of numbers of small order intervals in 
C. Up to factors of order one, we have in 2D and 4D, 
respectively: 



is-f^nC] =N-2Ni 
h 



and 



iV - A^i + 9iV2 - I6N3 + 8N4 , 



(13) 



(14) 



where N is the number of elements in C and Ni is the 
number of (i + 1) element inclusive order intervals in C. 

Because B is the most non-nonlocal of the operators 
in the family, the action S[C] is a sum of contributions 
each of which is not close to the value of the Ricci scalar 
at the corresponding point of the continuum approxima- 
tion. However, one might expect that if the curvature 
is slowly varying on some intermediate scale, which we 
might as well call Ik , the averaging involved in the sum- 
mation might perform the same role of suppressing the 
fluctuations as the smearing out of the operator itself so 
that the whole action S[C] is a good approximation to 
the continuum action when Ik is the appropriate size. 

There are many new avenues to explore. Can we use 
these results to define a quantum dynamics for causal 
sets? In 2D is there a relation with the Gauss-Bonnet 
theorem? Can we analytically continue the action in an 
appropriate way [l^l to enable Monte-Carlo simulations 
of the path sum? What sort of phenomenology might 
emerge from such actions? To answer this latter question, 
we need to know how big Ik must be so that the action 
S[C] is a good approximation to the Einstein-Hilbert ac- 
tion of the continuum SsHig]- In a rough estimate 
is reported that in dimension 4, Ik >> (PlY^^. Tak- 
ing L to be the Hubble scale, that would mean that in 
the continuum regime, only spacetimes whose curvature 
was constant over a scale {PLy^^ would be able to have 
an approximately local fundamental action. One might 
expect therefore that the phenomenological IR theory of 
gravity that could emerge from such a fundamental the- 
ory would be governed by an effective Lagrangian 
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where hi and 62 are of order 1, 60 is set to its observed 
value, and where Ik varies with epoch and today is much 
larger than the Planck scale. The phenomenological im- 
plications of these ideas remain to be explored. 



We end by pointing out that these results have a rele- 
vance beyond causal set theory as they provide a "proof 
of concept" for the mutual compatibility of Lorentz in- 
variance, fundamental spacetime discreteness, and ap- 
proximate locality. 
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